Universal Scaling of the Neel Temperature of Near-Quantum-Critical 
Quasi-Two-Dimensional Heisenberg Antiferromagnets 
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We use a quantum Monte Carlo method to calculate the Neel temperature In of weakly coupled 
S — 1/2 Heisenberg antiferromagnetic layers consisting of coupled ladders. This system can be 
tuned to different two-dimensional scaling regimes for T > Tn. In a single-layer mean-field theory, 
Xs D {Tn) = {z2 J') , where Xs D is the exact staggered susceptibility of an isolated layer, J' the inter- 
layer coupling, and Z2 = 2 the layer coordination number. With a renormalized 22, 22 — > fez, we find 
that this relationship applies not only in the renormalized-classical regime, as shown previously, but 
also in the quantum-critical regime and part of the quantum-disordered regime. The renormalization 
is nearly constant; hi ~ 0.65 — 0.70. We also study other universal scaling functions. 

PACS numbers: 75.10.Hk, 75.10.Jm, 75.30.Kz, 75.40.Mg 
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Antiferromagnets with effectively low dimensionality, 
consisting of weakly coupled chains (quasi-lD) or layers 
(quasi-2D), offer unique opportunities to study quantum 
mechanical collective behavior. A multitude of quasi-lD 
and quasi-2D antiferromagnetic compounds have been 
discovered, or deliberately designed, and they exhibit a 
wide range of ordered and disordred phases. At the same 
time, new and improved experimental techniques enable 
increasingly sophisticated studies of their properties. It 
is thus possible to test in detail microscopic quantum 
spin models and theoretical quantum many-body con- 
cepts, such as quantum-critical scaling 0. Motivation 
for studying these systems often come from phenomena 
directly associated with low dimensionality. Real materi- 
als, however, almost always have some three-dimensional 
(3D) couplings that cannot be completely ignored at low 
temperatures. The ways in which these couplings change 
the physics, e.g., leading to phase transitions or dimen- 
sional cross-overs [2| , are also governed by the physics of 
the ID or 2D units. Studies of 3D effects can therefore 
also provide important insights. 

In this Letter we investigate the finite-temperature 
Neel transition temperature Tn of a quasi-2D 5=1/2 
Heisenberg model consisting of layers of coupled lad- 
ders. In the absence of inter-layer couplings, the Mermin- 
Wagner theorem dictates that the system can have long- 
range order only at T = 0. The 2D Heisenberg model 
with spatially isotropic nearest-neighbor couplings J has 
an ordered ground state 0, 0- At low temperatures, 
in the renormalized classical (RC) regime, its spin cor- 
relation length is exponentially divergent [f|. Systems 
with a coupling pattern favoring formation of nearest- 
neighbor singlets e.g., coupled two- leg ladders 0, 
can be tuned through a quantum phase transition into a 
quantum disordered (QD) state. This T = transition 
and its associated T > quantum-critical (QC) scaling 
regime have been studied in detail, using field-theoretical 
approaches 01 and quantum Monte Carlo (QMC) simula- 
tions jdilSllldl- Some QMC studies of the Neel transition 
in systems of weakly coupled spatially isotropic layers 
have also been reported. Sengupta et al. studied Tn and 



the dimensional cross-over in the specific heat ^lj. Ya- 
suda et al. studied Tn more systematically at very small 
ratios a — J' /J of the inter- and intra-layer couplings, 
for S = 1/2 as well as higher spins [T^] • 

We here carry out QMC calculations of the coupled- 
ladder system to test three recently proposed scaling 
functions relating Tn and various 2D and 3D staggered 
susceptibilities 13, 14] . We focus in particular on inter- 
ladder couplings for which the system is near-quantum- 
critical in the absence of inter-layer couplings. Our main 
finding is that two of the scaling functions are almost 
constant and do not change appreciably between the RC 
and QC regimes. In particular, the coordination number 
renormalization introduced by Yasuda et al. changes from 
k 2 « 0.65 in the RC regime [H to fc 2 « 0.68 in the QC 
regime. This minute change implies that the Neel order- 
ing takes place almost exactly at the same temperature 
at which the planes start to correlate appreciably, and 
that these correlations are almost completely governed 
by the magnitude of the static staggered susceptibility of 
the planes. The nature of the fluctuations, RC or QC (for 
which the dynamic susceptibilities are completely differ- 
ent does not play a major role. 

In a single-layer mean-field theory (also referred to as 
RPA) |l5j , the 3D couplings of a layer I = are taken into 
account by a static staggered magnetic field arising from 
the ordered moments of the two adjacent planes I = ±1. 
The self-consistent Neel temperature is then obtained by 
solving the equation 



xl D {T N ) = (z 2 aY 



(1) 



where z = 2 is the layer coordination number and xi D 
is the staggered susceptibility of a single layer, the T de- 
pendence of which is known from studies of the quantum 
nonlinear a model. In the RC regime, T < Anrp s 



X ^(T)cxTe 4 ^/ T , 



(2) 



where p s is the spin stiffness. Using numerically exact 
QMC results for \ 2 S D an d T N , Yasuda et al. [3 found 
that the mean-field expression Q accurately captures 
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the aCl dependence of Tn, if z is replaced by a renor- 
malized coordination number z^k^ 

xl D (T N ) = (k 2 z 2 a)-\ (3) 

Moreover, the renormalization, k 2 ~ 0.65, was found to 
be independent on the spin S. This intriguing result 
prompted Hastings and Mudry to carry out a detailed 
renormalization group (RG) study of the anisotropic 
O(N) non- linear sigma model Instead of a constant 
coordination number renormalization, they argued that 
the quantity 

F 1 = (k 2 z 2 )- 1 = a X 2 s D (4) 

is a universal function of x = c[Xn<!; 2D (Tn)] _1 when a <C 
1. Here £ 2D is the correlation length of a single isolated 
layer and c the spinwave velocity. They concluded that 
the reason for the near constant k 2 is that the single 
layer is in the RC regime for all S at low T, whence x 
is exponentially small and F%(x — > 0) is constant. In the 
QC and QD regimes F\ should approach other constant 
values. A quantity involving the susceptibility x(Q) °f 
the full 3D system at wave-vector Q = (ir, ir, 0) was also 
introduced 

F 2 = ax(7r,7r,0). (5) 

To leading order in an 1 /N approximation, Praz et al. [hH ] 
found that F 2 — 1/4 in all regimes when a — » 0. This 
prediction should be easier to test experimentally because 
it involves only properties of the actual quasi-2D system. 
They also proposed a third universal quantity 

F 3 = aS(ir,Tr,0)T N \ (6) 

where 5(Q) is the static spin structure factor. This func- 
tion was shown to distinguish between the RC, QC, and 
QD regimes already at the N = oo level. 

Since the a — > values of F± , F 2 , and F3 were eval- 
uated at the N = go level or including only order- l/N 
corrections, significant higher-order corrections to these 
results were expected |l4j . Unbiased numerical results 
would therefore be useful. The previous QMC results by 
Yasuda et al. for k 2 = [2F(x — > 0)] _1 in the RC regime, 
k 2 rs 0.65 |l2j, falls between the N — 00 and l/N values; 
k 2 = 1/2 and 1.01 14]. F 2 and F 3 have not yet been 
calculated in the RC regime, and none of the predictions 
have been tested against numerical results in the QC and 
QD regimes. 

The universal constants could be very useful for ex- 
tracting inter-layer couplings experimentally This mo- 
tivates us to carry out large-scale QMC simulations of 
the coupled-ladder system, where the individual layers 
can be tuned through a quantum-critical point. We can 
then obtain numerical results for F1-F3 in all three 2D 
temperature regimes. 

The hamiltonian for coupled Heisenberg layers consist- 
ing of two-leg ladders is 

h = Ji J2 s i ■ s i + J 2 J2 Si • s j + J3 s * • s ^ ( ? ) 

<i,J>2 (i,i>3 




q c q=V J l 1 



FIG. 1: Finite-temperature regimes for a 2D layer of cou- 
pled ladders with inter-ladder coupling q. The curves indi- 
cate cross-over temperatures. We study T > properties of 
a quasi-2D system (coupled layers) with q at and close to the 
2D critical point q c and at the isotropic point q — 1. 



where denotes a pair of nearest-neighbor spins in 

the same ladder, (i,j} 2 in different ladders of the same 
layer, and in adjacent layers. We define the cou- 

pling ratios q — J 2 / J\ and a = J3/J1. 

The quantum phase transition of the single layer (a = 
0) has been studied by Matsumoto et al. ■ The critical 
coupling q c = 0.31407(5). The T > cross-overs for an 
isolated layer are shown schematically in Fig. ^ In the 
limit a — > 0, the quantum-critical coupling of the quasi- 
2D system approaches q c of the single layer. In our study 
we focus on values of q close to the 2D quantum-critical 
point, choosing q = 0.25, 0.30, 0.31407 = q c , and 0.33. 
We also consider the previously studied case q = 1 [l2T| . 
to calculate also F 2 and F3 deep inside the RC regime. 
We have obtained results for a in the range 10 -3 to 1. 

We use the stochastic series expansion (SSE) QMC 
method |l6j | to study periodic lattices with L x L y L z spins, 
with L x = L y = L up to 128. To take into account, at 
least partially, the fact that £ X} y 3> £ z when a <C 1, we 
use aspect ratios L/L z up to 16. To determine the Neel 
temperature, we use the finite-size scaling of the spin 
stiffness constants in the three different directions, 
fi = x,y, z, of the spatially anisotropic lattice. This ap- 
proach was previously taken in Ref. ,11]. For fixed as- 
pect ratio, the stiffness at Tjv should scale as We 
thus locate the point at which Lp^T) becomes asymp- 
totically size-independent (extrapolating crossing points 
for different size L to L —> 00). For q = 1, we use T)v 
from Ref. F° r the calculations of the 2D staggered 
susceptibility xi D {T) we have used L up to 800. 

The a dependence of the Neel temperature for the dif- 
ferent q values is shown on a log-log scale in Fig. El When 
q < q c , there is a minimum value a c of the inter-layer cou- 
pling below which the system cannot order — a c (q < q c ) 
is the line of 3D quantum-critical points and a c (q c ) = q c . 
In our results for q > q c we see a down-turn of Tn as a de- 
creases, reflecting the 3D quantum-critical point. From 
our limited 1ow-Tn data we can only roughly extract two 
points on the critical line; a c (q = 0.30) w 0.001 and 
a c (q = 0.25) w 0.006. 

For q > q c we can solve the mean-field equation 
with the RC form J2J of the correlation length, giving, to 
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FIG. 2: Neel temperature as a function of inter-layer coupling 
at different inter-ladder couplings q. The q = 1 results are 
from Ref. [T^. The solid (black) and dashed (red) lines shows 
the expected a — > forms in the RC and QC regimes. 




FIG. 3: Staggered susceptibility Xs vs - temperature for 2D 
systems composed of coupled two-leg ladders with different 
inter-ladder coupling ratios q. The dashed line shows the 
asymptotic QC power-law behavior, Eq. (1101 ) . 



leading order in a, 

Tjv(a) oc -[ln(a)]- 1 . (8) 

As shown in Fig. |5J for q = 1 this form does not yet 
apply at a = 1CP 3 , but it should be the correct form for 
a — ► 0. Yasuda et al. presented an empirical formula that 
works well also at higher a [l^ . For q = 0.33 we should 
also approach the RC form when a — > 0, but here we 
instead observe an almost perfect power law in the whole 
range of a > 10 -3 . However, the exponent is not the 
one expected in the QC scaling regime (discussed further 
below), and we expect the behavior to eventually cross 
over to the log form. 

In the QC regime, the 2D staggered susceptibility takes 
an asymptotic T — > power-law form, 

X 2 s D (T)^T- 2+ \ (9) 

where 77 ss 0.038 ^3] is t ne correlation function exponent 
of the 3D 0(3) Universality class. Using this form in the 
mean-field equation we get a corresponding power- 
law behavior of TV; 

T N (a) oca 1/(2 -"). (10) 

In Fig. |21 we can see that this quasi-2D quantum-critical 
form accurately describes the results for q = q c be- 
low a ~ 10 -2 . For larger a, Tn is still in the high- 
temperature regime where the behavior is influenced by 
non-universal lattice effects :5, 8J. The two q values close 
to q c , for which the reduced coupling \g_g c \/g c « 0.05, 
are already too far from the critical point to observe any 
distinct (asymptotic-form) QC behavior before the cross- 
overs occur. 

Following Ref. |12| , we study the coordination number 
renormalization 

k 2 (a) = pax^CTV)]" 1 = (2F 1 )-\ (11) 

In Fig. |3| we show our QMC results for the staggered sus- 
ceptibility of the isolated 2D layers. Using these results 



and the T N data shown in Fig.|2 we obtain the results for 
k 2 shown in the upper panel of Fig. El F or q = l, Yasuda 
et al. found k 2 « 0.65 for a < 0.1 [lj- We here show 
q = 1 results obtained with their listed Tn values and our 
own results for x 2D - The resulting k 2 agree with the pre- 
vious results. Surprisingly, we hardly see any change in 
k 2 when going to the near-critical systems, except some 
small differences when a > 0.1. At lower a, k 2 is only a 
few percent larger for q sa q c than at q = 1; k 2 (q c ) « 0.68. 
Even for our q < q c points, k 2 remains close to this value, 
even though Tn is seen crossing over into QD behavior 
in Fig. |21 For the lowest a considered for q = 0.25 and 
0.30 we see a slight increase in k 2 , but the effect is barely 
statistically significant. Note that for q < q Cl k 2 is not 
defined for a < a c (q). 

In the RG study by Praz et al. [3, different expres- 
sions for Fi (k 2 ) were obtained in saddle-point approxi- 
mations for the RC, QC, and QD regimes. No numerical 
values were given, however, except for k 2 — 1/2 in the 
RC regime. Corrections to the constant behavior was ex- 
pected (and calculated to order 1/N in the case of the 
RC regime, then giving k 2 — 1.01). Furthermore, sig- 
nificantly different constants were expected for the three 
regimes. The near constant k 2 w 0.65 — 0.70 we find 
here for such a wide range of q values, spanning all three 
temperature regimes, is thus quite remarkable. 

We now turn to the second scaling function, Eq. JSJ. 
As seen in the middle panel of Fig.0] we obtain an almost 
constant F 2 = ax(7r, 7r, 0) « 0.22 — 0.23 for all q and for a 
wide range of a. In the N = 00 approximation, F 2 = 1/4 
in all temperature regimes |l4j , remarkably close to what 
we find here. However, also in this case the actual values 
in the RC, QC, and QD regimes were expected to differ 
markedly once 1/N and higher corrections are included. 

The third scaling function, Eq. 10, distinguishes be- 
tween the RC. QC, and QD regimes already at the 
N = 00 level [II. In the RC regime S(Q) = T X (Q) 
to leading order Q and thus F3 = F 2 = 1/4. Our re- 
sults for the RC regimes (q = 1), shown in the bottom 
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FIG. 4: Inter-layer coupling dependences of the quantities Q, 
and J|)J| for different inter-ladder couplings q. The dashed 
lines show our extracted constant values for the QC regime. 

panel of Fig. 0] are slightly lower than the predicted value 
for a < 0.1. There is also still a decreasing trend as a 
decreases and we cannot reliably extract the asymptotic 
constant RC value (for a — 10~ 3 our calculations for 
q = 1 are not completely size converged and we therefore 
do not show them here). For q = 0.33, which also should 
give RC behavior for a — * 0, the results are still quite far 
from the q = 1 curve for all a, but the decreasing trend 
is consistent with the same asymptotic value. For F$ we 



also see clear differences in behavior in the three regimes. 
There are distinct cross-overs from QC to RC or QD be- 
havior. The results for q = q c ,0.30, and 0.33 all fall on 
the same universal QC curve for a down to 0.05, below 
which the q = 0.25 curve splits off. The q — q c and 0.3 
curves coincide to even lower a. The asymptotic value at 
q c is w 0.27. For q < q c we expect a divergence at a c (q), 
as 5(7r,7r,0) must converge to a constant when q < q c 
and Tn — > 0. We see clear signs of this divergence. 

In conclusion, we have presented results for the quan- 
tities ax^ D = 1/2&2, ctx(n,7T, 0), and aS(n,7r,0)/T, at 
T = Tn for a quasi-2D system of coupled ladder. For 
weak inter-layer coupling, a — > 0, these quantities have 
been predicted to take different universal constant values 
in the RC, QC, and QD regimes 0,0]. We have inves- 
tigated the dependence on a for 10~ 3 < a < 1, with the 
goal of extracting the constants and investigate the a > 
corrections. We find a remarkably stable value of the co- 
ordination number renormalization ki and cvx(7r, 7r, 0): 
For a < 0.1 they are almost independent on a and do 
not change appreciably between the RC, QC, and QD 
regimes. Significant differences in the three regimes were 
anticipated based on the previous RG study of the non- 
linear a model by Praz et al. |l4j|. Only in aS(ir, 7r, 0) 
do we see distinct differences. It would be useful to ex- 
tend the calculations to still lower inter-layer couplings, 
but reaching significantly below a = 10~ 3 with QMC 
requires prohibitively large lattices. 

The almost constant Fi and F2 imply that the corre- 
lations between layers is predominantly governed by the 
magnitude of the static staggered susceptibility of the 
layers. The range of temperatures for which the system 
is 3D critical is almost negligible, regardless of the nature 
of the 2D fluctuations — RC or QC — that initially lead to 
correlations between the layers. 
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